Mixed quantum state detection with inconclusive results 
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We consider the problem of designing an optimal quantum detector with a fixed rate of inconclusive 
results that maximizes the probability of correct detection, when distinguishing between a collection 
of mixed quantum states. We develop a sumcient condition for the scaled inverse measurement to 
maximize the probability of correct detection for the case in which the rate of inconclusive results 
exceeds a certain threshold. Using this condition we derive the optimal measurement for linearly 
independent pure-state sets, and for mixed-state sets with a broad class of symmetries. Specifically, 
we consider geometrically uniform (GU) state sets and compound geometrically uniform (CGU) 
state sets with generators that satisfy a certain constraint. 

We then show that the optimal measurements corresponding to GU and CGU state sets with 
arbitrary generators are also GU and CGU respectively, with generators that can be computed very 
efficiently in polynomial time within any desired accuracy by solving a semidefinite programming 
problem. 

PACS numbers: 03.67.Hk 



I. INTRODUCTION 

Quantum information theory refers to the distinctive 
information processing properties of quantum systems, 
which arise when information is stored in or retrieved 
from quantum states. A fundamental aspect of quan- 
tum information theory is that non-orthogonal quantum 
states cannot be pcrfcctly distinguished. Therefore, a 
central problem in quantum mechanics is to design mea- 
surements optimized to distinguish between a collection 
of non-orthogonal quantum states. 

We consider a quantum state ensemble consisting 
of m positive semidefinite Hermitian density operators 
{pi, 1 < i < Tri} on an n-dimensional complex Hilbcrt 
space 7í, with prior probabilities {pi > 0,1 < i < 
to}. For our measurement we consider general positive 
operator-valued measures [jï], 0], consisting of positive 
semidefinite Hermitian operators that form a resolution 
of the identity on Ti. 

Different approaches to distinguishing between the 
density operators Pi have emerged. In one approach, 
the measurement consists of to measurement operators 
which are designed to maximize the probability of correct 
detection. Necessary and sufficient conditions for an op- 
timum measurement maximizing the probability of cor- 
rect detection have been developed (| [|, || . Closed-form 
analytical expressions for the optimal measurement have 
been derived for several special cases f|, 0, |[ [H], • In 
particular, the optimal measurement for pure and mixed- 
state ensembles with broad symmetry properties, referred 
to as geometrically uniform (GU) and compound GU 
(CGU) state sets, are considered in pü| . Iterative pro- 
ccdurcs maximizing the probability of correct detection 
have also been developed for cases in which the optimal 



measurement cannot be found explicitly P, |Ï2|| . 

More recently, a different approach to the problem 
has emerged, which in some cases may be more use- 
ful. This approach, referred to as unambiguous quan- 
tum state discrimination, combines error free discrim- 
ination with a certain fraction of inconclusive results 
|| ü, || || |g , Ü H p. The bàsic idea, pioneered 
by Ivanovic Jl5|, is to design a measurement that with 
probability (3 returns an inconclusive result, but such that 
if the measurement returns an answer, then the answer is 
correct with probability 1 . In this case the measurement 
consists of to + 1 measurement operators corresponding 
to to + 1 outcomes, where to outcomes correspond to 
detection of each of the states and the additional out- 
come corresponds to an inconclusive result. Chefles (ïlfl 
showed that a necessary and sufficient condition for the 
existence of unambiguous measurements for distinguish- 
ing between a collection of quantum states is that the 
states are linearly independent pure states. The optimal 
measurement minimizing the probability (3 of an incon- 
clusive result when distinguishing between GU and CGU 
pure-state sets was considered in [JhJ , and was shown un- 
der certain conditions to be equal to the equal-probability 
measurement (EPM). 

An interesting alternative approach for distinguishing 
between a collection of quantum states, fixst considered 
by Chefles and Barnct [^ÏJ and Zhang et al. |2^] for pure- 
state ensembles, and then later extended by Fiuràsek and 
Jezek 
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to mixed-state ensembles, is to allow for a cer- 
tain probability of an inconclusive result, and then max- 
imize the probability of correct detection. Thus, in this 
approach, the measurement again consists of to + 1 mea- 
surement outcomes; however, now the outcomes do not 
necessarily correspond to perfect detection of each of the 
states. Indeed, if the quantum states are mixed states 
or linearly dependent pure states, then perfect detection 
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of each of the states is not possible 
by allowing for inconclusive results, a 



I. Nonetheless, 
righer probability 
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of correct detection can be obtained in comparison with 
the probability of correct detection attainable without 
inconclusive results. 

Necessary conditions as well as a set of sufHcient condi- 
tions on the optimal measurement operators maximizing 
the probability of correct detection subject to the con- 
straint that the probability of an inconclusive result is 
equal to a constant (3 were derived in |g3[ , using Lagrange 
multiplier theory. It was also pointed out in p3[ that ob- 
taining a closed form analytical solution to the optimal 
measurement operators directly form these conditions is 
a difficult problem. 

In this paper we extend the results of |2^] in several 
ways. First, using principies of duality in vector space 
optimization, in Scction fill we show that the conditions 



derived in 23 are both necessary and sufHcient. We also 
show that the Lagrange multipliers can be obtained by 
solving a reduced size semidefinite programming prob- 
lem. This approach lends itself to cfficient computational 
methods which are guaranteed to converge to the global 
optimum. 

Second, we derive a general condition in Section ^ un- 
der which the scaled inverse measurement (SIM) is opti- 
mal. This measurement consists of measurement opera- 
tors that are proportional to the reciprocal states associ- 
ated with the given state ensemble, and can be regarded 
as a generalization of the EPM to mixed-state ensembles. 

Third, we develop the optimal measurement for state 
sets with broad symmetry properties. Specifically, in Sec- 
tion |v| we consider GU state sets defined over a finite 
group of unitary matrices. We obtain a convenient char- 
acterization of the SIM and show that the SIM operators 
have the same symmetries as the original state set. We 
then show that for a pure GU state set and for vàlues of 
(3 exceeding a certain threshold, the SIM is optimal. For 
a mixed GU state set, under a certain constraint on the 
generator and for vàlues of (3 exceeding a threshold, the 
SIM is again shown to be optimal. For arbitrary vàlues 
of (3, the optimal measurement operators corresponding 
to a pure or mixed GU state set are shown to be GU 
with the same generating group, and can be computed 
very efficicntly in polynomial time. 

In Section |VÏ| we consider CGU state sets [E3 , in which 



the states are generated by a group of unitary matrices 
using múltiple generators. We obtain a convenient char- 
acterization of the SIM for CGU state sets, and show that 
the SIM vectors are themselves CGU. Under a certain 
condition on the generators and for vàlues of [3 exceed- 
ing a threshold, the SIM is shown to be optimal. Finally 
we show that for arbitrary CGU state sets and for ar- 
bitrary vàlues of [3, the optimal measurement operators 
are also CGU, and we propose an efhcient algorithm for 
computing the optimal generators. 

It is interesting to note that a closed form analytical 
expression exists for the optimal measurement when dis- 
tinguishing between GU and CGU (possibly mixed) state 
sets with generators that satisfy a certain constraint, un- 
der each of the three approaches outlined to quantum 



detection, where in the last approach we assume that (3 
exceeds a certain threshold. Furthermore, as shown in 
Jj], Q and in Scctions |v| and VI, the optimal measure- 
ment operators corresponding to GU and CGU state sets 
are also GU and CGU respectively, under each one of the 
three outlined optimality criteria. 

Before proceeding to the detailed development, we pro- 
vide in the next section a statement of our problem. 



II. PROBLEM FORMULATION 

Assume that a quantum channcl is prepared in a quan- 
tum state drawn from a collection of given states repre- 
sented by density operators {pi, 1 < i < to} on an n- 
dimensional complex Hilbert space Tí. We assume with- 
out loss of generality that the eigenvectors of pi , 1 < i < 
to, collectively spanj37| Tí so that m > n. Since pi is 
Hermitian and positive semidefinite, we can express pi 
as pi = 4>i4>* for some matrix ój, e.g., via the Cholesky 
or eigendecomposition of pi [£5| . We refer to <f>i as a fac- 
tor of pi . The choice of <pi is not unique; if <f>i is a factor 
of pi, then any matrix of the form = 4>íQí where Qi is 
an arbitrary matrix satisfying QíQ* = I, is also a factor 
of pi. 

To detect the state of the system a measurement is 
constructed comprising to + 1 measurement operators 
{Lli, < i < to} that satisfy 



IL > 0, < i < to; 

5> = /. 



(i) 



i=0 



Each of the operators Ilj, 1 < i < m correspond to de- 
tection of the corresponding states pi,l < i < m, and 
n corresponds to an inconclusive result. We seek the 
measurement operators IL that maximizc the probabil- 
ity of correct detection, subject to the constraint that the 
probability of an inconclusive result is equal to a constant 
13 <l. 

Given that the transmitted state is pj , the probability 
of correctly detecting the state using measurement oper- 
ators {IL,1 < i < to} is Tr(pjHj) and the probability 
of a detection error is J27Li i^j ^(Pj^-i)- Therefore, the 
probability of correct detection is given by 



Pd = y^Tr^IL), 



(2) 



where pi > is the prior probability of pi, with J2íPí = lj 
and the probability of a detection error is given by 



m m 



i=l j=lj=£i 



(3) 



The probability of an inconclusive result is 

m 

Pi = ^PiTrQ^IIo) = Tr(An ) = (3, (4) 
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where for brevity we denote 



such that 



i=l 



(5) 



Our problem is to find the measurement operators 
{lli, < i < m\ that maximize Pd of (||) subject to 
the constraints (|Tj) and (Q). 

Note that since Tr(pi) = 1 for all i, 



Pd+Pe + Pi= Y.Pi^iPi) = !• 



(6) 



When Pg = the states are distinguished unambiguously 
so that if outcome i is obtained for some 1 < i < m, then 
the state is with probability one. It was shown in ]Ï3]] 
that with P < 1 we can choose measurement operators 
such that Pe = if and only if the state ensemble is 
a linearly independent pure-state ensemble consisting of 
density operators pi of the form pi = \^>i){<j>i\ for a set of 
linearly independent vectors \<j>%)- If the vectors \4>i) are 
linearly dependent, or if the ensemble is a mixed-state 
ensemble, then Pe cannot be equal to 0. Nonetheless, 
we may seek the measurement operators that minimize 
Pe, or equivalently, maximize Pd, subject to Pj = (3 for 
some (3 < 1. 

Equipped with the Standard operations of addition 
and multiplication by real numbers, the space B of 
all Hermitian n x n matrices is an n 2 -dimensional real 
vector space. As noted in f23| , by choosing an ap- 
propriate basis for B, the problem of maximizing Pd 
subject to (Q) and can be put in the form of a 
Standard semidefinite programming problem, which is 
a convex optimization problem; for a detailed treat- 
ment of semidefinite programming problems see, e.g., 
p6| , |2Ï] , |28| , p9| . Recently, methods based on semidef- 
inite programming have been employed in a variety of 
different problems in quantum detection and quantum 
information |, Q IS 13 ü Ü- B y exploiting the 
many well known algorithms for solving semidefinite pro- 
grams [^9), e.g., interior point methods |Q |2^, ^8|, the 
optimal measurement can be computed very efficiently in 
polynomial time. 

The semidefinite programming formulation can also be 
used to derive necessary and sufficient conditions for op- 
timality, which we discuss in the next section. 



III. CONDITIONS FOR OPTIMALITY 

Using Lagrange multipliers, it was shown in [^3) that a 
set of measurement operators {11^, < i < m} maximizes 
Pd subject to Pj = /3 for a state set {pi,l < i < m} 
with prior probabilities {pi, 1 < i < m} if there exists an 
Hermitian X and a constant 5 satisfying 



X > PíPí, 1 < i < m; 
X > 5A, 



(7) 
(8) 



(X-p t p t )n t = 0, 1 < i < m; (9) 
(X-5A)fio = 0. (10) 

It was also shown that (||) and ( |Ï0| ) are necessary condi- 
tions for optimality. 

In Appendix ^ we use duality arg uments similar to 
those used in || to show that (|7j)-(|l·0|) are necessary and 
sufficient conditions for optimality, so that a set of mea- 
surement operators H maximizes Pd subject to Pf = (3 
if and only if there exists an Hermitian X and a con- 
stant 8 satisfying (ff|)-(|Ï0|). Furthermore, we show that 
X and 5 can be determined as the solution to the follow- 
ing semidefinite programming problem: 



min tr(X) 
xeB,Seiz 



5(3, 



(11) 



where 1Z denotes the reals, subject to 



X 
X 



Pi Pi 
SA. 



1 < i < m; 



(12) 



The problem of (11)-(|Ï2]) is referred to as the dual prob- 
lem. 

Note that the dual problem involves many fewer de- 
cision variables than the primal maximization problem. 
Specifically, in the dual problem we have n 2 + 1 real dc- 
cision variables while the primal problem has (m + l)n 2 
real decision variables. Therefore, it is advantageous to 
solve the dual problem and then use (JoJ) and ([!(]) to de- 
termine the optimal measurement operators, rather than 
solving the primal problem directly. 

The necessary conditions (0) and (^J) together imply 
that the rank of each optimal measurement operator 
is no larger than the rank of the corresponding den- 
sity operator; see ^j. In particular, if the quantum 
state ensemble is a pure-state ensemble consisting of 
(not necessarily independent) rank-one density operators 
Pi = \4>í){4>í\, then the optimal measurement is a pure- 
state measurement consisting of rank-one measurement 
operators Ui = \p,i){fMÍ 

As pointed out in |2^] , obtaining a closed-form analyt- 
ical expression for the optimal measurement operators 
directly from the necessary and sufficient conditions for 
optimality is a difficult problem. Since (^ÏJ) is a (convex) 
semidefinite programming |2^, |28|, ^9 ^problem, there are 
very efficicnt methods for solving ( l·l|) . In particular, 
the optimal matrix X and optimal sealar S minimizing 
Tr(A) — S/3 subject to ( |Ï2| ) can be computed in Matlab 
using the linear matrix inequality (LMI) Toolbox (see 
H [Ï4| for further details). Once we determine X, the 
optimal measurement operators Ui can be computed in 
a similar manner to that described in H. 

A suboptimal measurement that has been suggested 
as a detection measurement for unambiguous quantum 
state discrimination between linearly independent pure 
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quantum states, is the EPM (ï|, ||, |§, in which the 
measurement vectors are proportional to the reciprocal 
states associated with the states to be distinguished. A 
general condition under which the EPM is optimal for 
distinguishing between pure quantum states was derived 
in [|14j. It was also shown that for GU state sets and 
for CGU state sets with generators satisfying a certain 
constraint, the EPM is optimal. 

In the next section we consider a generalization of the 
EPM to mixed quantum states, which we refer to as the 
scaled inverse measurement (SIM). We then use the nec- 
essary and sufhcient conditions for optimality to derive a 
general condition under which the SIM is optimal. In Sec- 
tions |y| and VI we consider some special cases of mixed 
and pure-state sets for which the SIM is optimal, and 
derive explicit formulas for the optimal measurement op- 
erators. 



IV. THE SIM AND THE OPTIMAL 
MEASUREMENT 

The SIM corresponding to a set of density operators 
{pi = fa<j>*,l < i < m} with eigenvectors that collec- 
tively span Ti and prior probabilities {pí,1 < i < to} 
consists of the measurement operators {Ei = pip*,0 < 
i < m} where 



(13) 



for some 7 > and Eo = / — X)í=i/ í í/ í í· Here ^ is 
the matrix of (block) columns ipi = ^fplfa. Note that 
since the eigenvectors of the {pi} collectively span H, 
the columns of the {ipj} also together span Tí, so vpí'* is 
invertible. From (|Ï3|), 

m / m \ 

£ IHÚ = 7 2 A- J Yl Mi A- 1 = 7'A- 1 , (14) 

t=\ \i=\ ) 

so that 



(15) 



If follows from (|T|) that the SIM operators satisfy (|ï|) if 
and only if 7 2 < À„ where {Ai, 1 < i < n} denote the 
eigenvalues of A = í'í'* and A„ = min Ai. 

In the case in which the prior probabilities are all equal, 



m = 7 ($$* 



1 < i < m, 



(16) 



where $ is the matrix of (block) columns fa. 

Since the factors fa are not unique, the SIM factors 
Pi are also not unique. If pi are the SIM factors corre- 
sponding to fa, then the SIM factors corresponding to 
<f>'i = 4>iQi with QiQ* = I are p! i = p t Qi. Therefore, 
although the SIM factors are not unique, the SIM oper- 
ators Si = pip* are unique. 



The SIM corresponding to a pure-state ensemble \fa) 
consists of the measurement vectors \pi) — 7A _1 |^j), 
where \tpi) — \fpl\fa)- If m addition 7 = %/Anj then the 
SIM vectors are equal to the EPM vectors JÏJ] . 

The probability of an inconclusive result using the SIM 

is 

Pi = Tr(AEo) = Tr(A) - 7 2 Tr(J) = 1 - n 7 2 . (17) 
Therefore to satisfy (H), 



7 



1-/3 



(18) 



Since we must also have 7 2 < A„ we conclude that j3 
must satisfy 



A 



P > 1 - n\ n = /?, 



(19) 



For linearly independent pure quantum states it was 
shown in 1 14 that the SIM with 7 = \f\^ minimizes the 
probability of an inconclusive result subject to the con- 
straint that Pfj — 1 for state sets with strong symmetry 
properties. The smallest possible probability of an incon- 
clusive result in this case is fi = /3 m i n . It turns out that 
for a large class of state sets, including those discussed in 
|Ï4|| , the SIM also maximizes Pjj subject to Pj ~ (3 for 
P > /3 m in- From the necessary and sufhcient conditions 
for optimality discussed in Section III it follows that the 
SIM is optimal if and only if the measurement operators 
Hi = PiP* ,1 < i < m and Ho = Eo defined by ( [Ï3|) and 
@ satisfy (g) and @ for some Hermitian X and con- 
stant 5 satisfying (^) and (||). A sufhcient condition for 
optimality of the SIM is given in the following theorem, 
the proof of which is provided in the Appendix. 

Theorem 1. Let {pi = </>íç!>*,1 < i < m} denote 
a collection of quantum states with prior probabilities 
felí < * < m}- Let {Ei = pip*,0 < i < m} with 
{p t = 7 A" 1 V'i,l < i < m} and E = / - Sí™ 1 S í 
denote the scaled inverse measurement (SIM) operators 
corresponding to {ipi = ^/pï4>i,l < i < to}, where 
7 2 = (1 — (3)/n, A = <!"!'* and ^ is the matrix with block 
columns ipi . Let X n — min Ai where Ai are the eigenval- 
ues of A. Then the SIM maximizes Pd subject to Pi = (i 
for j3 > /3 m in with P m - ln = 1 — n\ n if for each 1 < i < m, 
(l/^)p*ipi = ip*^ 1 'ipi = ctl, where a is a constant in- 
dependent of i. 

It is interesting to note that the condition of Theo- 
rem |l| is identical to the condition given in Theorem 1 
of [H| for the least-squares measurement, or the square- 
root measurement, to maximize the probability of correct 
detection when Pj = 0. 

As we expect, the condition ip* A" 1 ^ = al does not 
depend on the choice of factor fa. Indeed, if fa = faQi is 
another factor of pi with Qi satisfying QiQ* = I, and if 
í'' is the matrix of block columns ^ = y/pïfai = ^/pïfaQi, 
then it is easy to see that (^-)* (*'*'*) _ Ví = al if and 
only if ^*A _1 -0i = al. 
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For a pure-state ensemble consisting of density opera- 
tors pi = \(f>i){(j)i\ for a set of vectors \4>i), (íAíIA -1 ^) 
is the ith diagonal element of P = 1 í/*A~ 1 ^f = 
^"(vl/vl/*)- 1 ^. The matrix P is just the orthogonal pro- 
jection onto JVfè)- 1 , where J\í(^>) is the null space of 
If the vectors \<pi) are linear ly independent, then 
Mi^) 1 - = {0} so that P = I and (V^A" 1 ^») = 1 for all 
i. It therefore follows from Theorem [ï] that for a pure- 
state ensemble consisting of linearly independent state 
vectors, the SIM maximizes Pjj subject to Pj = j3 for 
any [3 > /3 min . 

If the state pi — 4>i4>* is transmitted with prior prob- 
ability pi, then the probability of correctly detecting 
the state using measurement operators Sj = //í/í* is 
PiTr(p*(j)i(j)*iii) = Tr^*^^/-^)- It follows that if the 
condition for optimality of Theorem ^ is met, then the 
probability of correctly detecting each of the states pi 
using the SIM is the same. 

For a pure-state ensemble consisting of states \(f>i) with 
prior probabilities pi, the probability of correct detection 
of the ith state is given by \(pi\ïpi)\ 2 . Since {pi\i})i) = 
7(V , i|A -1 |i/> i ) > for any set of weighted vectors \ipi), 
(p,i\4>i) is constant for all i if and only if |(/ií|V>i)| ' s con ~ 
stant for all i. Therefore, we may intèrpret the condition 
in Theorem [ï] for pure-state ensembles as follows: The 
SIM is optimal for a set of states \<pi) with prior proba- 
bilities pi and for (3 > /3 m i n if the probability of detecting 
each one of the states using the SIM vectors is the same, 
regardless of the specific state chosen. 

In the remainder of the paper we use Theorem ^ to 
derive the optimal measurement for mixed and (not nec- 
essarily independent) pure-state sets with certain sym- 
metry properties. The symmetry properties we consider 
are quite general, and include many cases of practical 
interest. 



V. GEOMÈTRIC ALLY UNIFORM STATE SETS 

In this section we consider geometrically uniform (GU) 
pij state sets in which the density operators pi are de- 
fined over a group of unitary matrices and are generated 
by a single generating matrix. We first obtain a conve- 
nient characterization of the SIM for GU state sets, and 
show that under a certain constraint on the generator the 
SIM is optimal when (3 > /3 m i n - In particular, for (not 
necessarily independent) pure-state ensembles the SIM 
is optimal. We then show that for arbitrary GU state 
sets and arbitrary vàlues of (3, the optimal measurement 
is also GU, and we develop an efficient computational 
method for finding the optimal generators. 

Let G — {Ui,í < i < m} be a finite group of m unitary 
matrices Ui. That is, Q contains the identity matrix /; if 
G contains Ui, then it also contains its inverse Uf = U*; 
and the product UiUj of any two elements of Q is in Q 
ü- 

A state set generated by Q using a single generating 
operator p is a set S — {pi — UipU*, Ui S G}- The group 



G is the generating group of S. Such a state set has strong 
symmetry properties and is called GU. For consisteney 
with the symmetry of S, we will assume equiprobable 
prior probabilities on S. 

If the state set {pi, 1 < i < m} is GU, then we can 
always choose factors <pi of pi such that {(f>i — Ui<j), Ui G 
G} where cf> is a factor of p, so that the factors (f>i are also 
GU with generator </>. In the remainder of this section we 
explicitly assume that the factors are chosen to be GU. 



A. Optimality of the SIM for GU States 

For a GU state set with generating group G, com- 
mutes with each of the matrices U t G % Conse- 
quently, T = ($$*) _1 also commutes with Ui for all i, 
so that from ( fïq ) 

Pi = ÍT4>i = jTU i( f> = jU t T^ = U if i, 1 < i < m, (20) 

where 



p = 7($$* 



(21) 



It follows that the SIM factors pi are also GU with gener- 
ating group G and generator p given by (21). Therefore, 
to compute the SIM factors for a GU state set all we need 
is to compute the generator p. The remaining measure- 
ment factors are then obtained by applying the group G 
to p. 

From (p0|) we have that 



(1/7)/^ = 



i 



7V TO 



p*U*Ui(j> 



1 



7V TO 



p*4>, 



(22) 



where (f> and p are the generators of the state factors 
and the SIM factors, respectively. Thus, the probabil- 
ity of correct detection of each one of the states pi using 
the SIM is the same, regardless of the state transmitted. 
This then implies from Theorem ^ that for a (not neces- 
sarily independent) pure-state GU ensemble the SIM is 
optimal when (3 > /3 m in- For a mixed-state ensemble, if 
the generator cj) satisfies 



*($$*)- 



al 



(23) 



for some a, then from Theorem |ï| the SIM is again opti- 
mal. 



B. Optimal Measurement for Arbitrary GU States 



If the generator cf> does not satisfy (|23|), or if (3 < (3 m \ n , 
then the SIM is no longer guaranteed to be optimal. 
Nonetheless, as we now show, the optimal measurement 
operators that maximize Pd subject to Pi = (3 for any 
(3 are GU with generating group G ■ The corresponding 
generator can be computed very efficiently in polynomial 
timc. 
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Suppose that the optimal measurement operators that 
maximize 



J({U i }) = ^2Tv(p i U i ), 



(24) 



subject to 



Pi({TU}) = 1 - ^Tr ( jr pilli = 0, (25) 



are ±T and let J = J({IT}). Let r(j,i) be the mapping 
from X x X to X with X = {1, . . . , m}, defined by r(j, i) = 



fc if U*Ui = U k . Then the measurement operators IIV 



UjTl rUtl) U*, l<ï<mandn^=í- X)2=i n i for an y 
1 < j < m are also optimal. Indeed, since 11; > 0, 1 < 
K ra and £™ i H» < ^ n.p ) > 0, 1 < i < m and 

E = u s i J n, et? < 1/; = /. (26) 
t=i Ví=i / 

Using the fact that pi — UipU* for some generator p, 

m 

i=l 
m 

= E Tr (^ fc *n fc [/ fc ) 

m 

= E Tr tó) 

i=l 

= J. (27) 

Finally, 

(m \ I in 

E = Tr E UjUipUÍUjUr^, 

i,s — l J Vz,s— 1 

(m 
E ^Wn* 

/ m 

= Tr E | , (28) 



so that from @, Pi({fi\ j) }) = Pj({Üí}). 

Since the measurement operators ITV are optimal for 
any j, it follows immediately that the measurement op- 
erators {rii = (V m ) J2T=i i i — * — m l an( i no = 



/ — y^j—i li; are also optimal. Now, for any 1 < i < m, 

^ m 

ïï « = -E r <"< ; 

= -E^^w 

fe=i 

= Ui^jru* k íi k u^u* 



UTIU*, 



(29) 



where II = (1/m) XX=i U k U k U k . 

We therefore conclude that the optimal measurement 
operators can always be chosen to be GU with the same 
generating group Q as the original state set. Thus, to 
find the optimal measurement operators all we need is to 
find the optimal generator II. The remaining operators 
are obtained by applying the group Q to II. 

Since the optimal measurement operators satisfy Ilj = 
U.mi*, 1 < i < m and p t = U lP U* , Tr(p,IL.) = Tr(pH), 
so that the problem (|2j) reduces to the maximization 
problem 



maxTrípII), 



(30) 



where B is the set ofnxn Hermitian operators, subject 
to the constraints 

II > 0; 

m 

J2UiUU*<I; 



1 - Tr ( E UipUill = [3. 



(31) 



The problem of (|3C|) and (31) is a (convex) semidefmite 
programming problem, and therefore the optimal II can 
be computed very efhciently in polynomial timc within 
any desired aceuracy ^(| |2^, ^9|, for example using the 
LMI toolbox on Matlab. Note that the problem of @ 
and ( pl·| ) has n 2 real unknowns and 3 constraints, in con- 
trast with the original maximization problem (0) subject 
to (|ï|) and (^) which has mn 2 real unknowns and m + 2 
constraints. 

We summarize our results regarding GU state sets in 
the following theorem: 

Theorem 2 (GU state sets). Let S — {pi = 

UipU* ,Ui 6 G} be a geometrically uniform (GU) state 
set on an n-dimensional Hilbert space, generated by a fi- 
nite group Q of unitary matrices, where p — (fxf>* is an 
arbitrary generator, and let $ be the matrix of columns 
<f>i = Ui<f). Then the sealed inverse measurement (SIM) is 
given by the measurement operators = pip*,0 <i< 
m with 

Pi = Uip, 1 < i < m, 



7 



where 

M = 7 ($$*)-V, 

with 7 2 = (1 — 0)/n, and So = J — MiM*- ^ e STM 
/ias í/ie following properties: 

1. The measurement operators Ei, 1 < i < m are GU 
with generating group Q; 

2. The probability of correctly detecting each of the 
states pi using the SIM is the same; 

3. If = al for some a, then the SIM 
maximizes Pp subject to Pj = for [3 > 
1 — n\ n where X n is the smallest eigenvalue of 
(1/ra) Y^iL·i Pi > I n particular, if <f> = \<f>) is a vector 
so that the state set is a pure-state ensemble, then 
the SIM maximizes Pd subject to Pj = (3 for any 
[3 > 1-nAn. 

For an arbitrary generator <f> the optimal measurement 
operators IIj , 1 < i < m that maximize Pp subject to 
Pj = (3 for any [3 are also GU with generating group Q 
and generator H that maximizes Tr(pTl) subject to II > 
0,E£Li UiW* < I, and Tr(J2T=i UipUiïl) = 1-/3. 



VI. COMPOUND GEOMÈTRIC ALLY 
UNIFORM STATE SETS 

We now consider compound geometrically uniform 
(CGU) |24| state sets which consist of subsets that are 
GU. As we show, the SIM operators are also CGU so that 
they can be computed using a set of generators. Under 
a certain condition on the generators and for f3 > /3 m i n , 
we show that the optimal measurement associated with 
a CGU state set is equal to the SIM. For arbitrary CGU 
state sets and arbitrary vàlues of (3 we show that the opti- 
mal measurement operators are CGU, and we derive an 
efhcient computational method for finding the optimal 
generators. 

A CGU state set is defined as a set of density operators 
S = {pik = 4>ik4>* lk A < i < l, 1 < k < r } such that 
Pik = UiPkU* , where the matrices {Ui, 1 < i < 1} are 
unitary and form a group Q , and the operators 1 < 
k < r} are the generators. We assume equiprobable prior 
probabilities on S. 

If the state set {p lk , 1 < i < l, 1 < k < r} is CGU, 
then we can always choose factors <pik of pik such that 
{4>ik = Ui(f>k,í < i < 1} where <\> k is a factor of p k , 
so that the factors <pik are also CGU with generators 
{0fc, 1 _ k < r}. In the remainder of this section we 
explicitly assume that the factors are chosen to be CGU. 

A CGU state set is in general not GU. However, for 
every k, the matrices {</>ifc,l < i < 1} and the opera- 
tors {pikA < i < 1} are GU with generating group Q. 
Examples of CGU state sets are considered in . 



A. Optimality of the SIM for CGU State Sets 

With $ denotin g th e matrix of (block) columns 4>ik, 
it was shown in |7], [24] that $$*, and consequently T = 
, commutes with each of the matrices Ui G Q. 
Thus, the SIM operators are T^k = HikfJ-ik: 1 — i — h 1 < 
k < r with 

Pik = jT(f)ik = jTUiífrk = UiPk, (32) 

where 

Mfe =77> fe = 7(*$TVfc· ( 33 ) 

Therefore the SIM factors are also CGU with generating 
group Q and generators pk given by ( [33] ) . To compute the 
SIM factors all we need is to compute the generators pik- 
The remaining measurement factors are then obtained by 
applying the group Q to each of the generators. 
From @, 

P*k4>ik = p* k U*U l (j) k = p* k (j) k , (34) 
so that from Theorem |] the SIM is optimal if 

Vkfik = l4>t T( t>k =al, l<k<r, (35) 
for some constant a. 



B. CGU State Sets With GU Generators 

A special class of CGU state sets is CGU state sets with 
GU generators in which the generators {p k = (f) k (j>1,l < 
k < r} and the factors <f> k are themselves GU. Specifically, 
{4> k = V k (f>} for some generator <f>, where the matrices 
{Vk, 1 < k < r} are unitary, and form a group Q. 

Suppose that Ui and Vk commute up to a phase factor 
for all i and k so that U t V k = V k Uie je ^ where 9{i,k) 
is an arbitrary phase function that may depend on the 
indices i and k. In this case we say that Q and Q commute 
up to a phase factor and that the corresponding state set 
is CGU with commuting GU generators. (In the special 
case in which 9 = so that U\V k — V k Ui for all i, k, the 
resulting state set is GU @). Then for all i, k, <£<!>* 
commutes with UiVk 0, and the SIM factors pi k are 
given by 

Pik = jT4>ik = ~iTU l Vk(ï> = U l V k p, (36) 

where p — "/Tip. Thus even though the state set is not 
in general GU, the SIM factors can be computed using a 
single generator. 

Alternatively, we can express pi k as pi k — Uip k where 
the generators p k are given by 

Pk = V k p. (37) 

From ( |37| ) it follows that the generators p k are GU with 
generating group Q — {V k ,l < k < r} and generator p. 
Then for all k, 

Pl<f> k = P*V k *V k <f> = p*<\>. (38) 
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If in addition, 



Finally, 



H*4> = j<f>*T(f) = al 



(39) 



for some a, then combining (p4|), ( |3q ) and ( j39| ) with The- 
orem we conclude that the SIM is optimal. In partic- 
ular, for a pure-state ensemble, Jl*4> is a scalar so that 
([59|) is always satisfied. Therefore, for a pure CGU state 
set with commuting GU generators, the SIM maximizcs 
Pd subject to Pj = (3 for > f3 m i n - 



C. Optimal Measurement for Arbitrary CGU 
States 



If the generators (j> k do not satisfy (|35|), or if j3 < /3 m in, 
then the SIM is no longer guaranteed to be optimal. 
Nonetheless, as we now show, the optimal measurement 
operators that maximize Pd subject to Pj = are CGU 
with generating group Q. The corresponding genera- 
tors can be computed very emciently in polynomial time 
within any desired accuracy. 

Suppose that the optimal measurement operators that 
maximize 



J({Il ik }) =J2Y,Tr(p ik IL ik ), 



(40) 



i=i k=i 



subject to 



P/({H lfe }) = 1 - ^Tr í È ) = & ( 41 ) 

\ i,j=l k,s=l I 

are Tlik, and let J = J({Tli k }). Let r(j, i) be the mapping 
from X x X to X with X = {1, . . . , /}, defined by r(j, i) = s 

if UjUi = U s . Then the measurement operators Tl^jJ = 

Ujllry^kU* for any 1 < j ' < l are also optimal. Indeed, 

since U ik > 0,1 < i < 1,1 < k < r and J2\=i Z)fc=i Hifc - 

I, nife > 0, 1 < i < l, 1 < k < r and 

i=l fc=i fe=i / 

Using the fact that pi k — Uip k U* for some generators p k , 

l r 



i=l fe=l 
l r 

= ^^Tr(p fc f/;ÍT sfe C/ s 

s=l fe=l 
l r 

= y^y^Tr(;0i fc n ifc ) 

i=l fe=l 

= J. 



, ,s=l /c,t=l 

/ 



= Tr E E UjUiPkUÏUiU, 
yi,s=i fc,t=i 

(i r 
E E ^iPfc^n. 
i,s=l k,t=l 

= Tr í E É > ( 44 ) 

so that from ©, P/({ÍT^}) = P/({n ifc }). 

Since the measurement operators 11^ are optimal for 
any j, it follows immediately that the measurement op- 
erators {n ife = (í/OEí^n». 1 < * < í. 1 < k < r } 

and ILq = I — J^- fc are also optimal. Now, for any 
1 < i < 1,1 < k < r, 

- 1 - 

n lfc = y E ' -•' ' / 

1 - 

s=l 

= UiíÍ k Ut, (45) 

where Ü k = (í/l) £* s=1 í/ s *fl sfc C/ s . 

We therefore conclude that the optimal measurement 
operators can always be chosen to be CGU with the same 
generating group Q as the original state set. Thus, to find 
the optimal measurement operators all we need is to find 
the optimal generators {H k , 1 < k < r} . The remaining 
operators are obtained by applying the group Q to each 
of the generators. 

Since the optimal measurement operators satisfy Tli k = 
UiU k U* and p ik = U t p k U* , Tr(p lfc n ife ) = Tr(p fc n fc ), so 
that the problem (0) reduces to the maximization prob- 
lem 



™x VTr(p fc II fe ), 

llfctO — 



(46) 



k=l 



subject to the constraints 

n fe >, \<k<r- 



(43) 



L I 

i=l k=l 

\i=l k,l=l J 



(47) 
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Since this problem is a (convex) semidefinite program- 
ming problem, the optimal generators Tlk can be com- 
puted very efficiently in polynomial time within any de- 
sired accuracy |2^, g9), for example using the LMI 
toolbox on Matlab. Note that the problem of ( |46| ) and 
( ^7| ) has rn 2 real unknowns and r + 2 constraints, in con- 
trast with the original maximization (^|) subject to (|ï|) 
and (Q) which has Irn 2 real unknowns and Ir + 2 con- 
straints. 

We summarize our results regarding CGU state sets in 
the following thcorem: 

Theorem 3 (CGU state sets). Let S = {pik — 
UiPkU*,l < i < 1,1 < k < r} be a compound geo- 
metrically uniform (CGU) state set on an n-dimensional 
Hilbert space generated by a finite group Q of unitary 
matrices and generators {pk = (j)k<j)* k ,l < k < r}, and 
let $ be the matrix of columns <pik = Unpk- Then the 
scaled inverse measurement (SIM) is given by the mea- 
surement operators S^fc = pikp* ik , 1 < i < 1,1 < k < r 
and So = / - J2i k P^P* lk with 



where 



Pik = UiPk 



p k = 7(***) Vfe, 



and 7 2 = (1 — (3)/n. The SIM has the following proper- 
ties: 

1. The measurement operators Eifc, 1 < i < 1,1 < k < 
r are CGU with generating group Q; 

2. The probability of correctly detecting each of the 
states 4>ik for fixed k using the SIM is the same; 

3. If (<&<&* ) _1 0fe = al for some a and for 1 < k < 
r, then the SIM maximizes Pd subject to Pi = fi 
with /3 > 1 — n\ n where X n is the smallest eigen- 
value of (1/lr) J2i. k Pik- 

If in addition the generators {4>k — Vk4>, 1 < k < r} are 
geometrically uniform with UiVk = VkUie^^ 1 '^ for all 
i, k, then 

1- Pik — UiVkp. where p = 7(<i>$*) _1 so that the 
SIM operators are CGU with geometrically uniform 
generators; 

2. The probability of correctly detecting each of the 
states (f>ik using the SIM is the same; 

3. If cf)*^*)- 1 ^ = al for some a, then the SIM 
maximizes Pr> subject to Pj = /3 with fj > 1 — 
n\ n . In particular, if (f> = \(f>) is a vector so that 
the state set is a pure-state ensemble, then the SIM 
maximizes Pd subject to Pi = (3 with f3 > 1 — n\ n . 

For arbitrary CGU state sets the optimal measurement 
operators 11^,1 < i < 1,1 < k < r that max- 
imize Pd subject to Pi = f3 for any fj are CGU 



with generating group Q and generators Hk that max- 
imize X)fe=i Tr{pkHk) subject to Hk > 0, 1 < k < r, 

J2i, k UiU k U: < I, and IV (EL £*,i=i t ^* C W') = 
r(l-/3). 

VII. CONCLUSION 

In this paper we considered the optimal measurement 
operators that maximize the probability of correct detec- 
tion given a fixed probability (3 of an inconclusive result, 
when distinguishing between a collection of mixed quan- 
tum states. We first derived a set of necessary and sufh- 
cient conditions for optimality by exploiting principies of 
duality theory in vector space optimization. Using these 
conditions, we derived a general condition under which 
the SIM is optimal. We then considered state sets with 
a broad class of symmetry properties for which the SIM 
is optimal. Specifically, we showed that for GU state 
sets and for CGU state sets with generators that sat- 
isfy certain constraints and for vàlues of (3 exceeding a 
threshold, the SIM is optimal. We also showed that for 
arbitrary GU and CGU state sets and for arbitrary vàl- 
ues of (3, the optimal measurement operators have the 
same symmetries as the original state sets. Therefore, 
to compute the optimal measurement operators, we need 
only to compute the corresponding generators. As we 
showed, the generators can be computed very efficiently 
in polynomial time within any desired accuracy by solv- 
ing a semidefinite programming problem. 
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APPENDIX A: NECESSARY AND SUFFICIENT 
CONDITIONS FOR OPTIMALITY 

Denote by A the set of all ordered sets II = 
{nj^IU € B satisfying (0) and @j with (3 < 1, and 
define J(II) = Y^iLi Pí^ÍPí^U)- Then our problem is 



max JÍII). 

neA v ' 



(Al) 



We refer to this problem as the primal problem, and to 
any II € A as a primal feasible point. The optimal value 
of J(II) is denoted by J. 

To derive necessary and sufficient conditions for opti- 
mality, we now formulate a dual problem whose optimal 
value serves as a certificate for J. As described in ||, a 
general method for deriving a dual problem is to invoke 
the separating hyperplane theorem |36| , which states that 
two disjoint convex sets|39 can always be separated by a 
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hyperplane. We will take one convex set to be the point 
0, and then carefully construct another convex set that 
does not contain 0, and that captures the equality con- 
straints in the primal problem and the fact that for any 
primal feasible point, the value of the primal function 
is no larger than the optimal value. The dual variables 
will then emerge from the parameters of the separating 
hyperplane. 

In our problem we have two equality constraints, 
SZLo H« = I an d Tr(An ) = (3 and we know that 
J > J(II). Our constructed convex set will accord- 
ingly consist of matrices of the form — I+Y^iLo H where 
lli G B and ü; > 0, scalars of the form (3 - Tr(An ), 
and scalars of the form r — J(TÍ) where r > J. We thus 
consider the real vector space 

C = Bx1lxn = {(S,x,y): S£B,x,y£TZ}, 

where TZ denotes the reals, with inner product defined by 



((W, z, í), (S, x, y)) = Tt(WS) + zx + ty. (A2) 
We then define the subset Çl of C as points of the form 



i=0 



(A3) 

where Ui G B, U t > 0, r £ TZ and r > J. 

It is easily verified that SI is convex, and G" 
Í2. Therefore, by the separating hyperplane theorem, 
there exists a nonzero vector (Z, a,b) £ C such that 
({Z, a, b), (Q, c, d)) > for all (Q, c, d) £ Q, i.e., 



Tr U -/ 



+ E n JJ + H/3 - Tr (An )) 
+a Ir-f^PiTiiUiPi) j > 



(A4) 



for all Ui £ B and r £ TZ such that H > 0, r > J. 

As we now show, the hyperplane parameters (Z, a, b) 
have to satisfy certain constraints, which lead t o th e for- 
mulation of the dual problem. Specifically, (A4) with 
Ui = 0, r — ► J implies 



aJ > Tr(Z) - 6/3. 



(A5) 



Similar ly, ( |AJ ) with r = J + 1, üj = for j ^ i, = 
í|x)(x| for one value 1 < i < to where |x) G C n is fixed 
and t — > +oo yields (íc|Z — apiPi|a;) > 0. Since |x) and i 
are arbitrary, this implies 



Z > api/?; 



1 < i < TO. 



(A6) 



With r = J + 1, üj- = for j ^ 0, n = 
where |a;) G C™ is fixed and t — > +oo, (A4) yields 
— 6A|a;) > 0, which implies 



Z > bA. 



(A7) 



With U t = 0,0 < i < to, r -> +oo, (|AJ) implies et > 0. 
If o = 0, then @ yields Tr(Z) < b/3 < b and @ yields 
Tr(Z) > ò. Therefore we conclude that a > 0, and define 



X = Z/a, S = b/a. Then flA5[ ) implies that 
T(X,Í)< J, 



(A8) 



where T(X, 5) = Tr(X) - 5(3, (g6|) implies that AT > p.p; 
for 1 < i < to, and jA7| ) implies that X > SA. 

Let r be the set of X £ B, 5 £ TZ satisfying X > 
Pipi-, 1 < i < to and X > SA. Then for any X, 8 £ Y, 
U £ A, we have 



T(X, 5) - J(U) = Tr í U t {X - PlPl 
\i=i 

+Tr(n (A-<5A)) >0. 



(A9) 



Since X £ T, from (|A8|) and (|A9| ) we conclude that 
T(A, <5) = J. 

Th us w e have proven that the dual problem associated 
with ( |AÏ|) is 



min <r(A) — <5/3, 



(A10) 



subject to 



X > PíPí 
X > 5A. 



1 < i < to; 



(All) 



Furthcrmore, we have shown that there exists an optimal 
X, S £ r and an optimal value T = T(X, S) such that 

Let lli denote the optimal measurement operators. 
Then combining (|A9|) with T = J, we conclude that 



(A-^pOU = 0, 
(X-ÍA)fio = 0. 



1 < i < to; 



(A12) 



Once we find t he o ptimal X and 5 th at m inimize the 
dual problem ( AlOj ), the constraints (A12) are neces- 
sary and sufhcicnt conditions on the optimal measure- 
ment operators IT. We have already seen that these 
conditions are necessary. To show that they are sufh- 
cient, we note that if a set of feasible measurement opera- 
tors IT satisfies (|AÏ2|) , then J^Zi Tr ( n ;(^ - Pi Pi)) = 
and Tr ((X - SA)Ü ) = so that from <M), J{U) = 



T{X,6) = J. 



APPENDIX B: PROOF OF THEOREM |l| 

In this appendix we prové Theorem [ï]. Specifically, 
we show that for a set of states pi = <f>i4>* with prior 
probabilities P i, if (1 / 7) /i*^ = al, 1 < i < to, where 
Mi = 7( v I /l I / *) _ Vi = 7A _1 Vi are the SIM factors and 
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"0i = \fpï4>ii then there exists an Hermitian X and a 
constant S such that 

X > V'iV'*, 1 < í < m; (Bl) 

X > ÓA; (B2) 

(X - ipiíp*)^^* = 0, 1 < i < m; (B3) 

(X - 5A)(I -7 2 A-!) = 0. (B4) 



Let X = aA and 5 = a. Thcn (gg) and (|Bj) are 

immediately satisfied. Next, since al = ip%£í\fii = 
^* a -i/2 A -i/2^ í it follows that 



Multiplying both sides of (B5) by A 1 / 2 we have 



oA > Vi^*, 



(B6) 



which verifies that the conditions (Bl) are satisfied. 
Finally, 

{X - = c^AA" 1 ^ - «7^ = 0, (B7) 



al > A-^^A- 1 ^ 2 . 



(B5) so that the conditions (B3) are also satisfied. 
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